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1. Introduction 

Over the last decade, our theories of the early universe have been promoted from 
speculation to a field of intense scientific study. The most important developments 
in our knowledge concern the nature of the primordial curvature perturbation (jjj, which 
is believed to have seeded temperature variations in the cosmic microwave background 
(CMB). It is now understood that ( must have had a spectrum which was close to scale 
invariance on the scales probed by the CMB [21 [3l H]. 

Many proposals have been made to explain how a primordial perturbation with an 
almost scale-invariant spectrum could have been generated in the early universe. The 
most widely-studied candidate is an era of inflation that may have taken place at high 
energy O |6l [TJ [H [9l [TOl [H], where "inflation" is defined to be any epoch in which the 
scale factor a undergoes acceleration, d > 0. Under this condition, local regions of the 
universe are exponentially driven to spatial flatness, homogeneity and isotropy [12], and 
each light bosonic field acquires a perturbation generated by amplification of quantum- 
mechanical vacuum fluctuations [131 HH [13 US]- The spectrum of this perturbation 
is close to scale-invariance when the universe inflates at a rate a /a which is almost 
constant. The curvature perturbation observed in the CMB is supposed to be a model- 
dependent mix of these fluctuations, yielding anisotropics in the temperature of the 
microwave sky which are compatible with observation. Inflation apparently provides 
a natural framework in which one can simultaneously understand both the large-scale 
regularity of the universe and its small-scale irregularity. 

Inflation is not a single model, but rather a whole collection of scenarios which fit 
into the above framework. The only necessary ingredients are: (i) a specification of the 
field content, which allows a division into 'light' and 'heavy' fields; (ii) a background 
evolution a{t) which gives rise to d > with the Hubble parameter H = a/a slowly 
varying; and (iii) a rule for generating ( from the light bosonic fields. 

This prescription is rather general and implies that many models (perhaps 
with wildly different and mutually incompatible microphysics) may simultaneously be 
compatible with the observational data, since they may make equivalent predictions for 
the spectrum of (. Therefore, we must expect that it will be difficult to learn about the 
microscopic physics which was operative during the very early universe: it will almost 
certainly be insufficient simply to study the spectrum of (. In order to distinguish 
between wildly different models of the early universe it is necessary to find another 
source of information. 

Fortunately, any detailed model of the inflationary era does not merely predict the 
spectrum of it also implies a subtle but calculable network of correlations between 

f There are two primordial perturbations commonly encountered in the literature. The first of these 
is the comoving curvature perturbation, written TZ, which is proportional to the laplacian of the Ricci 
curvature of comoving spatial slices. On the other hand, the uniform density curvature perturbation C is 
proportional to the laplacian of the Ricci curvature on spatial slices of uniform density. On superhorizon 
scales, TZ and C are equivalent up to a convention for signs [1]. 
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the higher-order moments. These moments collectively measure the so-called non- 
gaussianity of ( and arise from interactions among the quanta of the (^-field and the 
other constituents of the early universe. Non-gaussian effects from interactions of ( 
quanta have been extensively investigated over the last few years, with the hope that 
observations of such effects may be able to discriminate between different models for 
physics in the early universe [I71[T8l[I3[20lini[22l|23l[2ai25l[26l|2^ 
ISllMlESllMlETllMllMlllQlllIlliailalll^ However, more is possible. Interactions 
do not only imply non-gaussian statistics in the three- and higher n-point correlation 
functions: they also imply quantum corrections to all correlation functions, and in 
particular the power spectrum or two-point function. It is possible that such corrections 
may be large in their own right, demanding that they be taken into account in accurate 
analyses of the observational data, as recently suggested by Sloth [46l |47] . Regardless 
of their exact magnitude, by searching for signatures of such quantum corrections in the 
power spectrum and correlating the results with predictions for non-gaussian statistics in 
the higher n-point functions we obtain a more sensitive test of physics during inflation. 

A second powerful motivation for studying loop corrections is a simple point of 
principle. The tree-level formula for the spectrum of ( is widely used to make predictions 
for the amplitude and scale-dependence of fluctuations generated in a very large class of 
early universe scenarios. Before deciding what degree of credence we should attach to 
any of these predictions, it is necessary to thoroughly investigate whether the tree-level 
amplitude is a genuine approximation to the full quantum result. 

In this paper, the prospects for detecting quantum corrections to the power 
spectrum are assessed in the simplest model of inflationary physics, that of inflation 
with a single scalar (p and arbitrary potential V{(j)). In view of the importance of 
accurate comparison with the precision measurements which are becoming available, this 
issue has already attracted considerable attention. Early work by Mukhanov, Abramo 
& Brandenberger |1HI US] and Abramo & Woodard [50] demonstrated that significant 
effects were possible (see also Unruh [51]). Later estimates of loop effects were made 
in a large number of models [5ai53l|5aiSllMll57llMl[59l[60l|6ll|6ai^ 
in EH EH EQl HH [12], even at two-loop order [73]. Recently, Sloth ^ ST] determined 
the full fourth-order action for Einstein gravity coupled to a scalar field and used this 
to estimate the one-loop correction to the power spectrum of scalar field fluctuations. 
Although one would naively expect the loop correction to be suppressed by a factor of 
(if/Mp)^ ~ 10^^°, where H is the Hubble parameter during inflation. Sloth's calculation 
yielded a significant cumulative effect (as large as 70% in some models) which could affect 
the precision determination of cosmological parameters from CMB experiments. The 
unexpectedly large size of the loop correction in this estimate is due to an amplification 
by A^, the total number of e-folds of inflation which occur. Since A^ can be very large in 
models where inflation begins at around the Planck scale, it may dramatically modify 
the predictions of naive dimensional analysis. 

This paper attempts to address these issues using a formalism similar to that applied 
by Sloth [inillT]. However, in contrast to previous analyses, the estimate is divided into 
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two parts. First, the one-loop correction to the power spectrum of field fiuctuations 
is computed soon after horizon exit, using the slow-roll approximation to control the 
calculation. This correction is not observable by itself; in a second step, it must be 
combined with other correlators of the fields to yield the one-loop correction to the 
power spectrum of the observable perturbation ( long after horizon exit. The correct 
combination can be computed using the 6N formula [TH [751 [El UHl [23]. This two-step 
process has several advantages. We shall see that the loop correction is generally afflicted 
by divergences at late times and on large scales. The 6N formalism naturally resums 
these late-time divergences into time evolution, which allows the slow-roll approximation 
to be kept under control. On the other hand, the divergences on large scales can be 
controlled by performing the calculation within a finite box. In analogy with the late- 
time divergences, it has recently been shown by Byrnes et al. [45] that these divergences 
can be resummed into spatial variation on large scales. 

The present paper is concerned with the technical issue of computing loop 
corrections go the power spectrum of field fluctuations. This calculation involves the 
application of standard methods from quantum field theory, adapted to the case of 
an expanding spacetime. On the other hand, the assembly of field correlators into ( 
correlators is an essentially classical calculation using the 6N formula. For clarity, this 
calculation will be presented separately elsewhere [77] . 

In ^ the background evolution and perturbation theory of the single scalar field 
are briefiy described. The perturbations are characterized (as in more complex cases) by 
cubic and higher self-interactions which involve the time derivative of the perturbation. 
This has important consequences for the calculation of quantum corrections. These 
corrections are introduced in ^ In §3.11 a path-integral expression for a general one- 
loop, single- vertex correction to the power spectrum is given in the Schwinger formalism, 
and in §3.21 the question of deriving a correct path integral expression for theories with 
time-derivative interactions is considered. In such cases the correct path-integral formula 
is well-known to contain a ghost field, whose quanta do not appear in physical states but 
circulate in the loops which give rise to quantum corrections. The Feynman rules for 
this theory are written down in §3.31 In §l]the assembled formalism is used to compute 
the leading radiative correction to the one-point function of the field. This is of interest 
in its own right, but also provides a simple setting in which some subtle features of 
the calculational machinery can be resolved. The one-loop correction to the two-point 
function is computed in ^ and a brief discussion is given in ^ 

^ is introductory and merely serves to fix notation. The reader who is mostly 
interested in the computation of the two-point function ((50(ki)50(k2)) may wish to 
skip directly to ^ and dispense with §^|3]-[11 These sections are largely dominated by 
the question of setting up a correct formalism in which the one-loop correction may be 
computed. 

Units are chosen throughout such that h = c = Mp = 1, where Mp^ = SnG is 
the reduced Planck mass. The metric convention is (— , +, +, +), and the unperturbed 
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background is written in cosmic time t as 

ds^ = -dt^ + a\t) dx2. (1) 

It is frequently more convenient to employ a conformal time variable, defined by 
Tj = J^dt' /a{t'). Indices labelling spacetime coordinates are chosen from the beginning 
of the Latin alphabet {a,b, . . .); indices labelling purely spatial coordinates are chosen 
from the middle of the alphabet (z, j, . . .). Where multi-field models are under discussion, 
the different species of light bosonic fields are labelled with Greek indices (a, /3, . . .). 



2. Inflation from a single scalar field 

2.1. The background evolution 

The simplest microphysical model capable of supporting an inflationary epoch consists 
of Einstein gravity coupled to a single scalar field with potential V{<p), which can be 
taken to be arbitrary except that it must allow inflation for some values of 0. The field 
(J) is known as the infiaton. The combined action for this system is 

S=-\j d'a;v^{i?-V'^0V,0-2V(0)}, (2) 

where g = det gab is the metric determinant, d^x = dt d^x is the product of the spacetime 
coordinate differentials, and R is the spacetime Ricci curvature. The background field (f) 
is taken to be spatially homogeneous and the background metric is parametrized by the 
scale factor a, given in ([1]). The evolution of a is determined by the Friedmann equation 

3H' = ^<P' + V{<P), (3) 

where H = d/a is the Hubble parameter, and obeys the homogeneous Klein-Gordon 
equation 

4> + 3H^ + V'{4>) = 0, (4) 

in which a prime ' denotes a derivative with respect to 0. The condition that inflation 
occurs is a > 0, or e < 1, where the parameter e is defined by 

H 

(6) 

Using Eqs. ([3])-(j4]) one can show that an equivalent definition is e = 0^/2if^. 

When e obeys the stronger condition e <C 1, the rate of change of is negligible 
in comparison with the expansion rate H. In this case one says that the field is 
slowly rolling. Although slow-roll is not mandatory for inflation to occur, the near 
scale-invariance of the power spectrum imprinted on the scales which are observed in 
microwave background experiments suggests that slow-roll was approximately satisfied 
if the CMB perturbation has an inflationary origin. When e <^ 1 applies, perturbation 
theory in e and related small quantities is known as the slow-roll approximation. In this 
paper, we compute all effects to leading order in e. 



One-loop corrections to a scalar field during inflation 



6 



2.2. Scalar perturbations 

Now consider small spatially-dependent perturbations in the infiaton, = (j)o + S(j){t,x.), 
where 0o is the homogeneous background evolution and 50 obeys the smallness condition 
I (50 1 ^ 1 00 1- Since dominates the energy density of the universe by assumption, any 
perturbation in will lead to a perturbation in the metric. These perturbations can be 
parametrized by a scalar (the lapse), a spatial vector A^* (the shift), and a spatial 
metric hij, 

ds^ = -N"^ dt^ + hij{dx' + N' dt){dx^ + N^dt). (6) 

Because of general coordinate invariance, not all choices of {N,N^, hij} lead to different 
configurations of the gravitational field. This redundancy is removed by fixing a gauge. 
We will choose to work in the spatially flat gauge, where h^j is given by its background 
value hij = a'^{t)6ij. Having done so, the metric functions and A^* are completely 
determined in terms of 50 by the constraints implicit in the Einstein equations. 
These can be obtained by inserting Eq. in the action ([2]). One obtains 

S=-l j dtd^x v^|Ar(V>V,0 + 2V)-i^(E^^E,,-E2 + 7r2)|, (7) 

where Eij = ^hij — V(jA^j) is the "momentum" associated with hij, Vj is the spatial 
covariant derivative compatible with hij, and tt = — N^Vjcj) is the field momentumlj] 
The equations of motion for the lapse and shift follow by varying 5" with respect to 
A^ and A^* respectively, and do not involve time derivatives. Therefore they are not 
evolution equations but constraints and can be solved algebraically: A^ and A^* are not 
propagating fields. Once A^ and A^* are known they may be substituted in to obtain 
a reduced action which depends only on 50. 
The A^ constraint is 

VVV.0 + 2V+ ^{E'^E,, - ^2 + vr^) = (8) 
and the A^* constraint is 

^^-^^})| = ^V,0. (9) 



One solves Eqs. ([H])-© order by order in 50. We write 

oo / oo \ oo 

AT = 1 + ^ a^, and AT, = V J ^ J + ^ Pmi (10) 

m=l \m=l / m=l 

where am, '&m and Pmi are all mth order in 50 and the Pmi are chosen to be divergenceless, 
so that diPmi for all m. The expressions necessary to compute S to third order in 50 

t We are adopting a convention, used throughout this paper, in which repeated spatial indices in 
complementary raised and lowered positions are contracted with the spatial metric hij , whereas a pair 
of repeated indices which both appear in the lowered position are contracted with the Euclidean metric 
5ij. Thus, a^bi — "^i j h^''a,ibj, whereas aibi — '^^aibi. Spacetime indices obey the usual Einstein 
convention, and always appear in complementary raised and lowered pairs which are contracted with 
the spacetime metric gab- 
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were given by Maldacena in the comoving slicing [12] and rewritten in the flat slicing 
for multiple fields in Ref. [21]. The expressions necessary to compute 5* to fourth order 
were obtained in the flat slicing by Sloth [iGj HT] in an approximation where all vector 
modes were absent, and given in generality in Ref. [78] . 

We work to leading order in the slow-roll approximation. At first order in 6(f) the 
leading terms are o(e^/^) , 



1 

2H 



2H 



and f3ii = 0. 



At second order in 6(j) the leading terms are o(e°) 



AH 



1 



(12) 



(13) 



where bracketed indices (■ ■ ■) are symmetrized with total weight unity and Urs is defined 
by 

T^rs = dr6(j)ds6(f) + 54>drds6(j), (14) 

with S = tr its trace in the Euclidean metric. Eqs. (ITTl)-( fT3l) can be inserted into 
the action, Eq. ([7]), after which one obtains an expansion of S in powers of 50. The first 
non-trivial term is quadratic. At o(e°) it is equal to 



^2 



- / dt d^x 
2 



-{d5ct)f 



There is a cubic interaction whose leading term enters at o(e^''^) 
written 



(15) 

[21j , which can be 



5, 



(16) 



The quartic term has leading terms of order o(e°) [78j. These terms correspond to 



^4 



4a 



AH 



The free field action 5*2 and the interactions {5*3, S4} as written here are all accompanied 
by terms of higher-order in slow-roll parameters, which we neglect. One must be careful 
to ensure that this approximation is accurate, and we will return to this at various 
points in the analysis (see also Ref. [77]). 



2.3. Expectation values 

The observables in this theory are expectation values of products of n factors of 
the perturbation 50, taken at a common time but at distinct spatial coordinates 
{xi, . . . ,x„}. It is often more convenient to work with momentum space expectation 
values which are obtained by taking Fourier transforms with respect to the Xj, giving 
k-space correlators which are functions of {ki, . . . , k„}. 
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At tree level the one-point expectation value vanishes, (^^(k)) = 0, since S(f> is 
by definition a perturbation in the comoving region under consideration. However, the 
gravitational background is time dependent since the scale factor a{t) varies with t, and 
therefore the vacuum state of the theory is changing continuously. This effect leads to 
gravitational production of inflaton particles [79] . Therefore we must expect a non-zero 
one-point function to be generated radiatively, reflecting the emergence of quanta 
from the vacuum. This issue is discussed in more detail in §1] below. 

The two-point expectation value deflnes the power spectrum P{k), 

(5(/.(ki)5(^(k2)), = (27r)35(ki + k2)P.{h). (18) 

The subscript denotes evaluation at the time when the k-mode under consideration 
left the horizon, which is unambiguous by virtue of momentum conservation which 
requires ki = ^2- At tree-level, P*{k) = H'l/2k^. It is often useful to work instead 
with the so-called dimensionless power spectrum, which is related to P{k) by the rule 
V{k) = k^P{k)/2TT\ 

The three-point expectation value defines the bispectrum, B(ki, k2, k^), and the 
four-point expectation value defines the trispectrum, T(ki, k2, ks, k4), 

(50(kO^0(k2)50(k3)) = i2nMJ2 ^^)B{kl, k2, k,), (19) 

i 

and 

(50(ki)50(k2)50(k3)<50(k4)) = {27rfSiJ2 k2, kg, k4). (20) 

i 

Typically, B and T are proportional at tree-level to and respectively, multiplied 
by a momentum-dependent form-factor [l8l|19l[20l[22l[80l[8ll[2ll[78l[82]. 



3. Quantum corrections 

3.1. Loop corrections from Schwinger integrals 

Schwinger's formalism. The appropriate formalism for computing expectation values 
in any quantum field theory was outlined by Schwinger [83]. Consider the vacuum 
expectation value of any observable O, observed at some time t^:, and computed in some 
theory with light scalar fields {0"}. By inserting a complete set of states at any time 
tjj > t*, this expectation value can be written 

{Q\0\Q), = j [d</.1(fi|0J^ = </^")(05' = </^10|fi)„ (21) 

where \VL) is the vacuum state at t — )■ — oo, the subscript indicates that the 
fields in the expectation value are evaluated at t^., and 0jj denotes evaluated at 
tjj. The integral J[d<y9"] is taken over all three-dimensional field configurations at tjj. 
Each factor in the product of transition amplitudes on the right-hand side of (12T!) 
can be expressed using the conventional Feynman path integral formula [HH [HHl [86] , 
= V9|0|fi) = /[d0"]^ O expiS*, where S is the action functional and the integral 
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is taken over all field configurations which begin in the state \Q) and end in the state 
= (p). These boundary conditions on (p are schematically denoted by the limits fl 
and (f attached to 



The interacting vacuum. In order to evaluate such integrals by the usual Feynman 
diagram expansion it is necessary to remove these boundary conditions, so that we 
integrate over all unrestrictedly. We follow the analysis of Weinberg |65] . To remove 
the restriction that the field must begin in the vacuum state one can integrate over all 
obeying an arbitrary boundary condition at t — )■ — oo, after multiplying the integrand 
by the vacuum wavefunctional, "^[ip] = (0(— oo) = tplQ). This has the desired effect of 
restricting the integral to field configurations which begin in the correct vacuum. The 
exact expression for \1/ depends on what we assume about but because the theory 
is supposed to be free as t — )■ — oo it must be a gaussian in the fields [HU |65]. Therefore 
we assume 



oc 



n 



exp 



5(q+r)^],(g)V^"(q)V^"(r) 



n 



exp 



for some set of weight functionals {^^(g)}, where [tpjQip) is a convenient abbreviation 
for the integral. The expectation value fl2Tl) can therefore be written [65] 

(n|0|(])oc (^ny[d^"]j|(^ny[d/]^jexp(i5[0_])nexp' 



(23) 



n/[d0X]n O exp(i5[0+])nexp 

where 'f denotes Hermitian conjugation, and the fields in the two path integrals 
have been differentiated by the addition of subscripts '+' and ' — '. The overall 
constant of proportionality is irrelevant. Since f l2^ requires an integral over final field 
configurations, it is possible to drop the restriction on the fields 0± at tj, provided we 
guarantee that the + and — fields for each species share a common value at this time. 
This can be accommodated by inserting a (5-function into the integrand which constrains 
the fields to agree [65] 



e 



1 2 



(24) 



where e is positive. Also, the action is real by assumption so the only effect of Hermitian 
conjugation in (l23l) is to flip the sign of the 15* term. 



Solution for propagators. Suppose that the action corresponds to a free field, so that it 
can be written S = (27r)~^ J d^/ci d^/c2 dti dt2 0(ti, ki)A0(t2, k2)/2 for some differential 
kernel A. Eq. f l23p can be written as an unrestricted path integral over the fields 
{0+, 0_} of the form 

.0-(t2,k2)J 



n /[d0^^dC])exp 



^^'^^^'^dt,dt2$:rr,(^i>kO 



(2vr)= 



12 



2P 



(25) 
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where we have assumed that there are no hnear couphngs among the various species, T 
denotes a transpose, and is the (2 x 2) kernel 



In Eq. ( 126|) . is the (5-function — t^) and 5joo is the (5-function 6(tj + oo). We will 
also occasionally use the notation 6ij = 6{ti — tj). The field propagator matrix for any 

particular species, consisting of propagators 6*-+, G } which connect the 

+ and — fields, is found by inverting the quadratic term given in fl25|) . 

dt2 d'k2 K,2 ( ) = i(27r)3<5i35(ki - ka) ( J ) • (27) 

The subscript '23' indicates that 6*++ is a function of times and momenta in the form 
G++(t2, ts, ks), etc.; and similarly for the other G. 

Eq. ( 127|) splits into coupled equations for G , G+_ and It will shortly 

become apparent that the doublets and {G , are to be regarded as 

forming complex conjugate pairs, so half of these equations are related to the other half 
by complex conjugation. In the application of interest, A is given to leading order in 
the slow-roll approximation by the laplacian of exact de Sitter space, 

Ai2 = j-j-{aH^2) + (ki ■ k2)a(5i2. (28) 
oti at 2 

Write = (27r)35(ki + ks)^^^^. The G++ equation reads 

gc- + 3i7«,)^G- H- - i,SMG?l = -^S.. (29) 

G+_ obeys the homogeneous version of (l29l) . whereas G obeys the complex conjugate 

of fl29l) and G_+ its homgeneous complex conjugate. In addition, Eq. (|27|) gives G+_ 
and G- the boundary conditions 

SiiG^+_ = 5ijGL^ and 5ijG1% = 5i«Gf+. (30) 

Homogeneous equation. Consider any solution, say G, to the homogeneous version 
of fl29l) . Any such solution is a function of the single variable ti, which after changing 
to conformal time i] can be written in the form G{r]) = Ck{v)/(^{v) some function 
({rj) to be determined. [The dependence of the mixed propagators on a second time 
argument, t2, enters only through the boundary conditions (!30l) .] The mode function (k 
must obey 

C + {k^l - 2u) - iaHf{2 - e)]a = (31) 

where k is the common magnitude of ki and k2, a prime ' denotes a derivative with 
respect to rj, and t satisfies 

' - » O- (32) 
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Eq. (PT]) is equivalent to the condition Q + {k"^ — {aHy {2 — e)}(k = almost everywhere, 
together with the boundary condition {Ck/o?) — ?■ as — )■ — oo. Heuristically, this 
boundary condition can be accommodated most naturally by redefining the range of r] 
to include some evolution in imaginary time, rj i— )■ 77(1 + u). Although strictly speaking 
this contour is singular, owing to the presence of the (5-function, it can be approached 
as a limit of regular contours. It will be seen below that when integrations over 77 are 
required the integrands in question are holomorphic. Therefore, any one of these regular 
contours suffices for calculation and we may as well take 77 h-> 77(1 + ii) for fixed l. This 
prescription was used in Refs. [HI [211 HH] to compute tree-level correlation functions of 
the {50"} in the interacting vacuum. 

Propagator matrix. One can now construct an explicit solution for which satisfies 
(in conformal time with arguments rji and 772) 

G-,(k„k2) = (2vr)35(k, + k2) X I , (33) 

where denotes complex conjugation, k is the common magnitude of ki and k2, and 
Ck{'ni,V2) is defined by 

Uvi,v2) ^ i ^^/^^;y\'' a(r]i)a(r/2), (34) 

a{7]i)a{r]2) 

in which W{f,g) is the Wronskian W{f,g) = fg' - gf. Note that i[W{QXk)]~^ is 
real and time- independent, in virtue of Abel's identity, but may depend on k. The 

propagator G is obtained by complex conjugation of Eq. (l33ll : the mixed propagator 

G+_ is obtained from a homogeneous equation and therefore is smooth at rji = ri2- The 
boundary condition fl30|) implies that it must satisfy 

G^2_(ki,k2) = {2n)'6{k^ + k2)Uvum) (35) 

and G_+ is given by its complex conjugate. [Note that there is no ambiguity in deciding 
which propagator should be assigned to a mixed pair (50+(50_), because the mode ( 
is always assigned to the argument of the -|- field, and (* is always assigned to the 
argument of the — field.] 

The above analysis was carried out for a single field, but where more than one 
species of light field is present similar results apply, with a mode function for each 
species which obeys a vacuum boundary condition of the form {(^/a"^) — )■ in the far 
past. The propagators which connect two fields of different species a and /3 then obey 
analogues of fl33|) and (l35l) with the function C,k replaced by a matrix If the fields 
do not couple linearly, then it follows that = 5"^^^. 

One-vertex, one-loop amplitudes. In the remainder of this paper, we shall be concerned 
with computing expectation values in which a set of n external fields, 0(k„), observed 
at some time r/* and carrying momenta {k„}, are paired with a single {n + 2)-valent 
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internal vertex with coupling constant g. Applying Schwinger's formula shows that the 
term in such an expectation value of leading order in g is given by 

i(2vr) J (^^)3(n+2) 5(}_^qjy^^dr7(?M, (36) 

where M is defined by 

M ^ (0^(ki) ■ ■ ■ 0^(k„)0|(qi) ■ ■ ■ 0t(qn)/+(qn+i)0+(qn+2)) 

- (0^(kO ■ ■ ■ 0^(k„)0l(qi) ■ ■ ■0i(qj0^_(qn+i)r (qn+2)) • (37) 

Greek indices label the species of fields, which are here allowed to run over field 
derivatives as well as the fields themselves; the issue of obtaining a correct path 
integral for theories with derivative interactions causes no difficulties for the purposes 
of Eqs. fl36|) - fl37|) . but will be taken up again in more detail in the next section. Any 
amplitude of the type given in ( l36l) -( l37l) is automatically of one-loop order, because the 
two field operators left over after all n external fields have been paired with n of the 
vertex fields must contract amongst themselves, leaving a single unconstrained integral 
over momentum. 

The time integral in fl36|) has been carried to some arbitrary late time rj^ which 
satisfies rj^ > rj^,. Using Eqs. (I33l)-(l35l) together with their complex conjugates in 
Eq. (136|) . it follows that the expectation value can be written 

i{2nf6{k, + ... + K) J ^ |_%r/C(^,^*)---ea,r7.)er(^'^) 

- i(2vr)35(ki + ■ ■ ■ + k^ J ^ £jve,riv.v*)---Civ,v*)eriv,v) 

+ permutations, (38) 

in which the second term is the complex conjugate of the first, and all permutations 
likewise assemble into complex conjugate pairs. Observe that the internal term ^q{r],r)) 
in the first line comes from pairing two + fields, whereas in the second line it comes 
from pairing two — fields. Eq. f lM|) shows that for the S(f) propagator, ^q{ri, rj) is real, so 
that ^q{ri,ri) and C,q{ri,ri) are in fact equal. 

The part of the integration over times between rj^, and t]^ has cancelled out, since 
in this region Eq. fl35p is given by the same expression as Eq. fl33p . whereas in the 
region rj < rj^ it is given by its complex conjugate. Note that for interactions which 
contain more than a single vertex and a single loop the process of deriving expressions 
such as (1551) using the path integral technology described above becomes increasingly 
cumbersome. For such interactions, some form of the diagrammatic operator formalism 
recently elaborated by Musso [87] is likely to prove superior (see also Ref. [55]). 

3.2. Theories with derivative interactions 

An important feature of the interactions (fT6|) and (ITTl) is that they include time 
derivatives of the perturbation, S<p [SI]. This means that the lagrangian can not be 
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written in the canonical form L{6(j), 6(f)) = ^6(f)A6<p + V{S(f)) (where the operator A is 
field independent), in which there is only a quadratic dependence on 6<p. As a result the 
textbook construction of the path integral formula based on L does not work. 

In the standard construction one identifies a momentum, vr, canonically conjugate 
to 6(j) and writes the lagrangian as a Legendre transformation of the hamiltonian function 

L{6(f),6^) =7t6^- H{6(j),7r). (39) 

In the quantum theory 6(f) and vr cannot be specified simultaneously. Since it is H 
that generates time evolution, when one constructs the path integral one naturally 
arrives at a functional integration that involves independent integrals over 50 and vr. 
If L depends at most quadratically on 50 then H depends at most quadratically on 
TT and the momentum integral can be performed immediately. This has the effect of 
setting the value of vr equal to the one stipulated by Hamilton's equations and results 
in the standard lagrangian path integral formula [H]. However, when H has a more 
complicated dependence on vr the momentum integral must be treated more carefully. 

The properties of lagrangians with derivative interactions have been studied 
extensively in the context of the non-linear a-model. (See, eg., Coleman [89]; a path 
integral treatment is given in Ref . , whereas the canonical approach was followed in 
Ref. [90].) After inspection of Eqs. (fT6|) -( !T7|l it is clear that no term contains as many 
as four time derivatives, although there are terms containing one, two or three. Let us 
parametrize a general action for a field 6 with arbitrary interactions containing as many 
as three time derivatives in form 

s = (271)' j dr^ {^-io.pre^ - iL^dO'^do^ - vi9) + + ^Jo.M^'^e^e'^'^ . m 

In order to keep this and subsequent expressions manageable, Eq. (140|) has been written 
in an abbreviated "de Witt" notation, where contraction over indicies implies not only a 
summation over species, but also an integration over momentum variables with measure 
d^A;/(27r)'^. With these conventions the object 5aj3 = 5(kQ, + k/3) is a "pseudo-metric" on 
k-space which is numerically identical to its index-raised counterpart, 5"'^ 13 Note that 
we have taken any interactions involving exactly two factors of 9 to be included with the 
kinetic term in '~^ap- Moreover, without loss of generality 7^/3 and ujap-y can be supposed 
to be symmetric under exchange of their indices. We assume that 7q,/3 is invertible, with 
inverse 7"^^. 

The momentum conjugate to 0" is tt^, 

TT, = = {27if U^pO^ + A, + u^p^e^eA (41) 

where we have used the assumed symmetry under index exchange to simplify this 
expression. In order to apply this formalism to the cubic and quartic interactions ( |T6l) 

f Although it is tempting to regard Sap as an object for raising and lowering indices, it is not a true 
metric because with our conventions, the object obtained by mixing indices, SapS^'^ , is not the identity 
operator {2TT)'^5{ka — k^), although it is proportional to it. 
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and (IT7|) it is only necessary to compute to 0(^^), where we formally assume that vr ~ ^ 
in order of magnitude. Since there are no three-derivative interactions in 5*3, this implies 
that u = 0{6) and it will be sufficient to work to leading order in u. To this order, the 
hamiltonian can be written 

^ = ^tAs^"''^"^/^ + k'2Tr)'L(sde''d9^ + ^(27r) V^A^A^j + (27r)V 
Z [Zn Z Z 

- ^^a;„;3,7"^^'^7""vr,7r,7r. - T^t^^t^x. (42) 
This hamiltonian can be used to construct a path integral for 9, giving 

j [d^° dTT^j] exp |i j dri (nj'' - ^) | ■ (43) 

The fields 6°" and vr^ are now variables of integration, and therefore independent, so we 
are free to redefine the momentum field by a shift, 

TTa ^ {2nf{TTa + Xa) (44) 

with Xa chosen to eliminate the term in Eq. fj43|) which is linear in tTq, 

x^ = ^^0e^ + x^ + u^^^e^e\ (45) 

This shift leaves the path integral measure [dvT/j] invariant. Having done so, one may 
rearrange terms to find a simplified path integral expression 

y"[drd7r^] exp|i(^e + S'gh)}, (46) 

where Se is the original 6 action f l40|) with all derivative interactions in their original 
form, and S'gh is an effective action for the "ghost" field vr, 

^gh = (27r)3 j (^-^^<-^n^Trp + uj^^.^^^'^^^t^^^r + ^c^./^-yT^^^'^T^'vr.vr^Tr,^ (47) 

The quanta associated with tt do not appear in physical states, although they couple 
to 9 and so affect its expectation values when loop corrections are taken into account. 
This explains why it has been permissible to ignore such ghosts in previous tree-level 
calculations; the vr integral makes no contribution to tree-level expectation values. Note 
that unlike the more familiar Fadeev-Popov ghost, the tt field is a spacetime scalar and 
not a spin- 1/2 fermion. 

Eq. ( 1471) is not yet in a form suitable for perturbative calculations. In particular the 
ghost kinetic term involves the inverse 7"''. This will be a complicated object even for 
relatively simple choices of 7^/3, but it is pointless to compute beyond 0(6*^) since the 
action to which we wish to apply this formalism [namely Eqs. ( [T5|) - (fT7|) ] was truncated 
at this level. For a canonically normalised scalar field in an almost-de Sitter spacetime 
7a^ can be written 

7„/3 = a25(k„ + k/3) + 2ri(k„, k^) + 2r2(k„, k^) + ■ • • (48) 

where the F^ are taken to be o(^™'), the factors of two have been inserted for future 
convenience, and '■ ■ ■' denotes higher order terms which have been omitted. 
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The identity operator with our conventions is 62, = {2'7i)^S(ka — k^). Therefore, the 
inverse 7"^ can be written 

= (27r)^ (^'^(ka + + ^i(k„, kp) + V^2(k„, kp)^ , (49) 

where the ipm are taken to be o{6"^). One can verify that the normahzation in Eq. (l49ll 
is correct, since when (l48l) and (l49l) are contracted together the 0(6'°) term becomes 

/ ^ + k/3) ■ (27r)6l5(k^ + k^) = {2n)'6{K - k^) = S^. (50) 

The o(^) equation imphes that ipi satisfies 

2 

?/'i(k<„ k^) = -— ri(-k„, -k^), (51) 



(52) 



whereas the o(^^) equation imphes that ip2 satisfies 

2 4 f 

^2(k„,k;3) = -— r2(-k^, -k^) - — / d^gri(-kc„q)ri(-q, -k/j). 

The ghost action can therefore be written 

+ ^^c^a/37'^""'^^"'^""7i-p^<x7r.| (53) 

The first term is o(^)^ and can be taken as the free-field part, whereas the remainder 
is 0(^'^) and can be taken as the interaction term. In this form, the ghost action is 
suitable for perturbative evaluation. 



3.3. Feynman rules for the interacting scalar/ghost theory 

We are now in a position to write down the Feynman rules for the 6(j) theory, including 
the effect of the ghost field. In this section we will not be obliged to carry out any of the 
complicated manipulations which characterized §3.2[ and so we will revert to a notation 
in which momentum labels and integrals are written explicitly. 

The propagators for the pure S(j) theory were written down in §3.1[ The free part 
of the ghost action can be inverted immediately to find the ghost propagator. For the 
+ fields this gives 

(7r+(?7i, ki)7T+{T]2, ks)) = -j^a{rii)a{ri2)S{T]i - r/2)5(ki + ks), (54) 

from which the propagator can be obtained by complex conjugation. Eq. ( I54|) is 

the propagator for a so-called static ultra-local field [HI [85]. Its k and rj dependence 
is constrained by the appearance of 5-functions, so the ghost does not propagate: its 
purpose is to provide corrections to the vertices of the theory which account for the 
presence of coincident time derivatives there. At one-loop order, we do not require the 
mixed propagator; the ghost field only appears in loops, but at one-loop the only mixed 
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contractions involve pairings of external fields with internal fields and these cannot occur 
for the ghost. 

In the one- loop, single- field vertex formula, Eq. ( l38l) . the ghost propagator always 
appears in the role of the propagator evaluated at equal arguments, associated with the 
internal momentum q. Note that the special simplification, which occurred for the (50 
propagator, where ^q{ri,ri) and ^g{r],r]) were equal, does not apply for the ghost field. 

It remains to identify the interaction terms V, A, Fi, r2 and u. Reading these off 
from Eqs. ([I6])-([I7D we obtaiifl 

^ = "'sf / ^^^w^ {n^*(%)} • ^) (55) 

r.(q,.q,) = ~a^^ I ll ntn^^.) (l + 2^(f|S5l) (57) 

r,(q„q,) = -.^/ !!^«(X:q.){n«^(q,)} 

^ f z(qi,q3) • z(q2,q4) , 3 or(qi,q3)or(q2,q4) , ^(uM^h^^a) 

I 4 2 7 2 2 '~ 4 

-(qi,q2,q3) = / || ^(E q.)^'^(q4)^^, (59) 

where = q^ + q^. The functions a and z are defined by 

(T(a, b) = a ■ b + 6^ (60) 

and 

z(a, b) = cr(a, b)a — cr(b, a)b. (61) 

These are the momentum-space counterparts of Eqs. (fT3l) - (fT^ . Note that as written, 
Eqs. (l57I) - (!58l) for F and Eq. (159|1 for u are not symmetric under exchange of their 
arguments. For F this is immaterial, because (HQ]) and (15^ show that it always appears 
in a symmetric contraction. On the other hand, u does appear once in an asymmetric 
contraction, namely Uaji-yO'^^^'^S'^'^T^u'^T- To avoid an unnecessarily tripling of the length 
of (!59|) we leave it in asymmetric form, carrying out an explicit symmetrization when 
computing amplitudes involving the asymmetric vertex. 

X In vl-2 of the arXiv version of this paper, and the version which subsequently appeared in JCAP, 
the interaction r2 contained a sign error in its third term, which propagated through the remainder 
of the calculation. This caused the coefficient of In k which appeared in Eq. (|105p to be incorrect. I 
would like to thank P. Adshead, R. Easther and E. Lim for bringing this error to my attention. This 
error was corrected in Ref. |92j : however, owing to a typographical error in that reference, the final 
coefficient of In k was again given incorrectly (P. Adshead, personal communication) . 
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Figure 1. Pure Scf) vertices. A dot on a scalar line entering a vertex shows that a time 
derivative is applied to the field at the point of interaction. In terms of Eq. (|40|) . the 
diagrams correspond to the vertices produced by (a) the potential V: (b) the A vertex; 
(c) the Fi vertex; (d) the T2 vertex; and (e) the uj vertex. 
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Figure 2. Scalar/ghost vertices. Solid lines represent the scalar field S4>, whereas 
dashed lines represent the ghost. A dot on a scalar line entering a vertex shows that 
a time derivative is applied to the field at the point of interaction. Time derivatives 
are never applied to ghost fields. In terms of Eq. (1531) the diagrams correspond to the 
vertices produced by (a) the ^pl^^'^ interaction; (6) the ^p2'^'^ interaction; (c) the ujOtttt 
interaction; and (d) the ujtt^ interaction. 



Diagrammatic representation. Eqs. ( |55l) - (!59l) lead to a rather complicated diagram- 
matic formalism in which the vertices produce a number of related terms, depending 
on the number of derivatives which apply to the lines entering the vertex. In order to 
keep track of these related contributions it is useful to introduce a refinement of the 
Feynman rules in which the lines of scalar propagators to which derivatives are applied 
are decorated with a dot. 

For the pure 50 vertices, the resulting diagrams are depicted in Fig. [H For the 
mixed (50/ghost vertices, the resulting diagrams are shown in Fig. [21 

4. The one-point function 

In §2.3[ we observed that at tree-level the one point function of 5(j) is zero, ((50(k)) = 0. 
This is not merely a question of convention; if the one-point function was not zero 
then so-called 'tadpole' diagrams such as Fig. ([3]) would mean that 6(j) quanta would 
emerge from the vacuum. Conservation of momentum forces such particles to condense 
in the zero-momentum mode, and the accumulation of such particles causes the classical 
background field to change. Such an instability implies that any perturbation theory 
based on the original unstable vacuum state would not give meaningful answers. This 
problem can be avoided by ensuring that the vacuum which we take as the basis of our 
perturbation theory is stable, at least at tree-level. 

In an inflationary universe, the emergence of 6(f) quanta from the vacuum is exploited 
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Figure 3. Instability of tlie vacuum due to condensation. 6(f) particles emerge from 
the vacuum (represented by the hatched condensate) in a zero momentum state. The 
accumulation of such particles changes the homogeneous classical field configuration 
associated with the vacuum. 



to produce small density fluctuations on superhorizon scales. Therefore we may expect 
to encounter some symptoms of vacuum instability when quantum corrections are taken 
into account. These symptoms manifest themselves as a radiatively generated one-point 
function, 

{Sm) = i2nf5ik)0, (62) 

where O 7^ is a dimensionless quantity. Although it is not the observable in which 
we are principally interested, the present section is devoted to a calculation of O. This 
is important for two reasons. The first is that it provides a consistency check on 5N 
calculations |2ll ESI El, IH [82l HU |45] which typically assume = 0, even beyond 
tree level. The second is that it allows us to develop some aspects of the calculational 
formalism in a simpler setting than the computation of the two-point function. 

4.I. Ghost diagrams 

Consider the one-point function associated with some wavenumber k. We aim to 
compute this at the time rj^ when k crosses the horizon, which is roughly defined by the 
condition —hrj^, = 1. Eventually conservation of momentum will force us to set k = 0, 
but in order to regularize the calculation we compute for finite k and then study the 
limit — 0. 

We deal first with the diagrams which contain a ghost loop. There is only one such 
diagram, which arises from the ipiir"^ coupling, 

{Sm)* ^ O- (63) 

This diagram makes a contribution to (50(k))* equal to 

Mk)). 2 -(2.)3«(k) £d, / ^(0)(l - it.)e-^ A (^1 + 2'-t^^ 

+ complex conjugate, (64) 

where the symbol indicates that (50(k)) contains the indicated contribution (among 
others), and in deference to the vacuum prescription outlined in §3.11 we should deform 
the contour of the rj integral to include some evolution in imaginary time for large {rjl. 
In this region the exponential factor is strongly decaying (cf. the discussion in §3.ip . 
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so there is very little contribution to the integral from very early times; if 77* is not too 
late, the integral receives its dominant contributions from times around horizon crossing, 
where rj ~ —1/k. We may therefore approximate the slowly varying factors H^,H and 
4>/H by their values at the time of horizon crossing, which are equal to and 
respectively. In this simple example the rj integral and the integral over the internal 
momentum q factorize, leaving a final result 

+ complex conjugate (65) 

where P^:{k) is the tree-level power spectrum evaluated at rj^. The object 5(0) is the rj 
delta-function evaluated at zero argument, and is badly divergent. In the present case, 
however, this is not material. The rj integral can be rotated to imaginary time, leaving 
a result which is purely imaginary. Hence, although divergent, this diagram makes no 
contribution to the one-point function. 

4-2. Pure 6(f) diagrams 

Now consider the pure 5(j) diagrams. 

The vacuum prescription and renormalization. In these diagrams, as above, early times 
make almost no contribution to the rj integral, so that slowly varying quantities such as 
H and (p/H can be evaluated at ?7*. A generic pure-50 contribution to O will then take 
the form 

0* P* f d^g 

O ^ iP*(/c) — — / di] - — -- e'^'^S, +complex conjugate, (66) 
4:H^ J_^ J {2ny 

where E is a k- and q-dependent quantity which is to be calculated. In evaluating E 
we will encounter instances where a 6(j) propagator begins and ends at the same vertex, 
giving it coincident time arguments. We will choose to set such a propagator equal to 

(50(qi,r/)50(q2,r/)) = {27rf 6 {q, + ci,)^{l + q'r^'), (67) 

where the exponential factors e^'^'^e""^^ have cancelled among themselves, and q is the 
common magnitude of q^ and q2. 

The discussion of vacuum boundary conditions in §3.11 emphasized that the fields 
which participate in the Schwinger path integral must be chosen to begin in the 
appropriate interacting vacuum, and that this could be achieved heuristically by 
deforming the contour of integration to include some evolution in imaginary time. In 
view of this, one may question whether ( l67l) is the correct choice, or whether it should 
be modified to reacjjj 

tt2 

(50(qi,r/)50(q2,r/)) = {27i)^{q^ + q2)77^|l - igr^l^e'^'^^'^^''), where Im(r/) > 0. (68) 



2q 



f It makes no difference if we allow to develop a small imaginary component in the prefactor |1 — igTyp, 
since the exponential term is so strongly decaying for large \ri\. 
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This would apparently have the very desirable effect of decoupling our prediction for 
{S(j){k)) from any details of the deep ultraviolet regime where g — )■ oo, because fl68|) is 
strongly decaying in this limit owing to the exponential factor. Therefore one might 
have some reservations that the choice of ( 167|) would introduce unphysical divergences, 
arising from an incorrect treatment of the vacuum. On the other hand, Eq. ( I68l) has 
the undesirable feature that it leads to a non-holomorphic integrand. This means that 
it would be necessary to rescind the possibility of contour rotation in evaluating the 
rj integral. A loop amplitude computed using ( l68l) would therefore depend sensitively 
upon the supposedly arbitrary value we assign to liai{ri), which in turn depends on the 
structure of the vacuum at past infinity. It is not clear how the resulting amplitude 
should be interpreted. 

This situation can be understood as follows. It was explained in §3.1l that the trick of 
contour rotation can be expected to account reliably for the vacuum boundary conditions 
only when the integrand is holomorphic. This is automatically true at tree-level, as 
argued in §3.11 and gives completely unambiguous results which are independent of the 
details of physics in the deep ultraviolet. At one-loop level the situation is different. If we 
adopt Eq. ( 1671) . then the loop amplitude depends on the ultraviolet parameter lm{ri). If 
we adopt the holomorphic expression Eq. ( 167|) for the propagator, then the result is free 
of any dependence on Im(?7) but there is instead a prospect of sensitivity to the details of 
ultraviolet physics from the q integral. Although we can choose from which source this 
sensitivity arises, we cannot evade it altogether — as we should expect. However, any 
divergences associated with the limit g — ?• oo can be subtracted by conventional methods 
of renormalization, after which the contour-rotation prescription gives a finite, contour- 
independent result which correctly incorporates the vacuum boundary conditions. In 
what follows we adopt this prescription. 

This leaves open the question of how the g-integral should be regulated. Since the 
Einstein action is supposed to be an effective theory of gravity for energies less than 
the Planck scale Mp ^ 10^^ GeV, and inflation is usually supposed to occur at energies 
at least a few orders of magnitude less than Mp, one might imagine applying a cutoff 
on the loop momenta of order the Planck scale. However this in itself is ambiguous 
since the Planck scale, unlike the speed of light, is not a Lorentz invariant and varies 
between locally inertial frames. In particular, the comoving Planck scale at a given 
instant rj is given by a{ri)Mp. Therefore a momentum cutoff of this form entangles 
the rj and q integrations, and for this reason it seems preferable to use a method of 
regularization, such as dimensional regularization, which does not depend on the explicit 
use of a cutoff. In the present paper we will compute expectation values using a fixed 
momentum cutoff both in the infrared and ultraviolet. In the case of the one-point 
function it may be checked that when the ultraviolet region has been discarded both 
dimensional regularization and a fixed momentum cutoff yield comparable predictions 
for the leading infrared divergences. 

This is sufficient for the purposes of the present paper, since it is the behaviour in 
the infrared rather than the ultraviolet which is of principal interest in a cosmological 
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context. Divergences in the ultraviolet come from the behaviour of the fields at high 
energies and small scales. Such small scale modes exist far inside the horizon, where the 
equivalence principle suggests that flat spacetime quantum field theory is expected to 
be a good approximation. The subtraction of these modes has recently been considered 
by Finelli et ai, who argue that no special treatment is required for the power spectrum 
[93] (see also Ref. [M])- On the other hand, the infrared behaviour comes from low 
energies and large scales, where the field modes are well outside the horizon. On such 
scales, flat spacetime field theory is a very poor approximation and we are obliged to 
take account of the gravitational background. 

This does not preclude the appearance of new ultraviolet divergences in our 
expectation values. Indeed, many of the integrals we shall encounter do contain 
ultraviolet divergences, of which the ultraviolet divergent quantity 5(0) which appears 
in Eq. ( l65l) is an example. Such divergences do not interfere with our ability to perform 
accelerator or laboratory particle physics experiments on earth, which are characterized 
by time- and length-scales that are small compared to the expansion time-scale and 
horizon length-scale of the universe. On such small scales the ultraviolet divergences 
we shall encounter (none of which are present in the 50 theory in Minkowski space) are 
presumably subdominant with respect to divergences from the pure matter theory and 
therefore do not interfere with our ability to perform terrestrial experiments, or with 
the success of the principle of equivalence. 



Zero-derivative interactions. We now return to the 50 diagrams. It is simplest to 
classify these diagrams according to the number of derivatives applied to propagators 
entering the vertex. 

There is a zero-derivative interaction from the vertex in Fig. (II]) (a), 

(50(k)). D O- (69) 

This term makes a contribution to S which equals 

S ^ - f ^ - ^) (1 - ikv){l + qW). (70) 

The term involving k ■ q is not rotationally invariant and disappears in the integral over 
q. Let us introduce a fixed ultraviolet cutoff A and infrared cutoff /i. Evaluating the q 
and ?7 integrals as described above gives 

Two derivatives, both derivatives on internal leg. The next class of diagrams contain 
two derivative operators, and divide naturally into two sorts: those where the derivatives 
are applied to both ends of the internal loop, and those where one derivative is applied 
to the loop but the other is applied to the external leg. 
The first sort give rise to diagrams of the form 
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{Sm)* ^ O; (72) 

such diagrams contribute an amount to S corresponding to 

ES-|(1-..,)(1 + 2::(H1ZM). (73) 
Evaluating the integrals by the method described above, one arrives at 



Sg--(l-igr^)(2 + 2 "^ ^ ) (76) 



a 3 - + ^" + /^')- (74) 

Two derivatives, single derivative on internal leg. The final class of diagrams contain a 
single derivative on the internal line, and apply the remaining derivative to the external 
leg. These diagrams are of the form 

mk))* ^ — o (75) 

and contribute to S according to the rule 

k'^ , , / cr(— k, q) 

- 1-igr^ 2 + 2^'^^ 
2q V T 

This class of diagrams makes a contribution to O* which equals 
4-3. Infrared behaviour 

We now collect terms from Eqs. (ITTIl . (17i|l and (1771) . The result shows only power law 
divergences, which could have been anticipated from the outset because the tadpole 
loop must be independent of k. Accordingly, the — t- limit is perfectly regular, and 
yields 

This is purely divergent. Whatever renormalization scheme we choose, all these power- 
law divergences must be subtracted. The result may leave a fc-independent remainder, 
but even if this is true the zero-momentum limit /c — )■ will be finite. 

In practice, infiation does not last for an indefinite number of e- folds and the region 
of the universe described by the infiationary patch will not be unboundedly large. One 
should identify the k = mode with the spatial average of 6(f) within this patch. If 
((50(k) is non-zero after renormalization, this spatial average can be absorbed into a 
redefinition of the background field 0(t) by enforcing the condition ((50(k)) = 0, as 
discussed (for example) in Refs. [621 [631 SHI HZ] • It follows that when 6(f) is defined in 
this way one may take O = 0, as usually assumed. 
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5. The two-point function 

We now turn to the central purpose of this paper, the computation of the leading loop 
correction to the two-point function ((50(ki)50(k2))*. 

It is first necessary to decide which classes of diagrams are to be included in the 
computation. In general, the one-loop correction to the two-point function of the 50 
will be given by a sum of diagrams, the leading terms of which are of the form 



loop correction D 




Eqs. f|T6|) - f|T7|) show that the leading contribution from the first diagram is 0(e°) in 
slow-roll, whereas the leading contribution from the second diagram is 0(e). Therefore, 
provided e <^ 1 and there are no large logarithms which can compensate for small slow- 
roll parameters, the expectation value will be dominated by the lowest-order part of the 
first diagram. This is opposite to the case considered by Sloth jlHl ST] , where a large 
logarithm was used to compensate for the smallness of e. In this regime the first diagram 
will be dominated by its subleading slow-roll part, and for consistency one should also 
take into account sub-leading slow-roll terms from the second diagram, and possibly 
from other sources. In the present paper, we wish to use the slow-roll approximation 
to simplify the calculation and therefore we will retain only the contribution from the 
leading part of the first diagram. The question of when this is a good approximation, 
together with a more general analysis of any possible large logarithms, will be postponed 
to another publication [77] . 



5.1. Ghost diagrams 

The relevant vertices here come from the ip2'n"n' and uj6(j)TnT terms in the ghost action. 
There is no contribution from the wtttttt interaction because this involves three ghost 
fields, which must appear in loops, and at one-loop order there is always one ghost field 
which is left unpaired. Therefore this term can be disregarded, although it would play 
a role in a two- or higher-loop calculation. To determine the ip2 contribution explicitly, 
consider Eq. (!52l) which gives ■02 in terms of the known functions Fi and We are 
computing to leading order in slow-roll, so the term involving can be discarded, 
because Eq. fl57|) shows that it is proportional to the slow-roll parameter e ~ 
whereas the leading terms in the fourth-order interaction are o(e°). 
The relevant ghost diagrams are 

(50(ki)50(k2))* 3 + . (79) 

Both these diagrams are purely imaginary and cancel between the -|--|- and 

propagators in exactly the same manner described in Ig] for the computation of the 
one-point function. 

The ghost diagrams have therefore entirely cancelled out in both the one- and 
two-point functions. This leads to expressions which agree with those reported in 
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Refs. jSni HT] . However, one should not immediately conclude that the ghost diagrams 
always sum to zero. Although this issue deserves more detailed attention, at two-loop 
order and above one can presumably expect the factors of i to combine to give non- 
vanishing contributions. This will apparently occur whenever there are an even number 
of ghost propagators in the diagram. 



5.2. Pure 6(f) diagrams 

As in the one-point calculation, it is convenient to classify the pure 50 diagrams 
according to the number of derivatives they contain. 



Single derivative. There are no zero-derivative interactions, because the gravitational 
interactions responsible for generating the vertices in Fig. [1] make no contribution to 
the potential at o{6(j)^), and the cubic contribution which is generated would make a 
contribution to the loop correction which is subleading in slow-roll. 

The first non-trivial term contains a single derivative, which can be applied to an 
internal or external line, 

(<50(ki)50(k2)). D . + . (80) 

As in the case of the one-point function, it is useful to parametrize the contribution 
each diagram makes to (50(ki)50(k2))* in terms of a function 11, which is defined by 

((50(ki)50(k2)), = i(27r)3p,(/c)2 J dr] J -0^ e'^'^^'H + complex conjugate, (81) 

where P^{kY is the square of the tree-level power spectrum, and the quantity 11 (to 
be calculated in this section) depends on the external momenta {ki,k2} and the loop 
momentum q. 

The class of diagrams where the derivative is applied to the external leg makes a 
contribution to 11 which corresponds to 



n.(l-..,)(l + ,V)^^(q-k.)^^ + 3^ 

+ [ki ^ k2], (82) 

where [ki ^ k2] denotes the same term with ki and k2 interchanged. The ratio 
cr(— ki, k2)/A;i2 is obviously singular when the momentum conservation condition ki -|- 
k2 = is enforced and must be treated carefully to avoid an unphysical divergence. 
Consider the non-singular quotient (T(a, b)/|a + bp where a and b approach k and — k 
respectively, 

a(k + ^,-k + 6) _ \.-{6 + e) + ^ + 6-e 

^TTp " efeo 5^ + e^ + 25-e 

This is not symmetric between 6 and e, because a is not a symmetric function of its 
arguments; as a result, the limits do not commute. Moreover, as 6 and e approach zero 
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the numerator of f l83p vanishes hnearly, as fast as 0{e,S), whereas the denominator is 
vanishing quadratically, hke 0(6^,5^). Therefore fl82|) is naively divergent. In fact, the 
value of ( 183|) depends on what is assumed about k ■ d and k ■ e; if we demand that the 
limit is approached along a sequence of vectors of magnitude = |k| then it follows that 
|k + ^| = |— k + e| = fc and therefore 

2 £2 

k ■ e = - and k d = . (84) 

2 2 ^ ' 

With this choice, Eq. ( l83l) evaluates to 1/2 and the limits become commuting. This 

prescription was used implicitly in §4.2 of Ref. [78], but there does not seem to be 

any compelling reason to demand that the limit is approached along such a specific 

sequence of vectors. Fortunately a catastrophic divergence is averted, since Eq. fl82|) 

requires symmetrization over the exchange ki ^ k2. The problematic term k • (<5 + e) 

is antisymmetric under this exchange and cancels out of the expectation value ( 182|) . 

leaving a finite limit. The result of this procedure gives the same answer as if we had 

adopted Eq. (^^, which can be regarded as a justification for the analysis presented in 

Ref. m- 

After performing the symmetrization over ki and k2 and integrating over q and t], 
this class of diagrams make a contribution to the two-point function of the form 



((50(ki)50(k2)), D (27r)35(ki + ka)-^ - — P In A; - — P + ■ ■ ■ (85) 



' Pin A; P 

7r2 V 16 120 

where '■ ■ ■' denotes ultraviolet power law divergences which have been omitted. 

Now consider the diagrams in which the derivative is applied to the internal line. 
Such diagrams contribute to 11 according to 

rr^n -7 ^2/1 • x /ki ■ q (T(-q, -kg) /c^ a(q, -q) 

^ : U D [1 - ikr]) [1 - iqri) i — ■ 7^ + ^-, ^ 

\ 4g |k2 + q|^ 8g |q — q| 

+ [ki ^ ka], (86) 

This class of diagrams contains a similar ill-defined ratio, cr(q, — q)/|q — q|. Consider 
Eq. (!83|) again, with k replaced by q. Although there is no longer any injunction to 
symmetrize over q ^ — q, the non-rotationally-invariant part q ■ (<5 + e) will vanish 
underneath the integral and does not give rise to any divergence. In order to assign a 
definite value to the remaining limit, we must assume something about S and e. Since 
q is merely a variable of integration and can be freely replaced by — q, we assume that 
o"(q, — q) is to be regularized by taking its symmetric part. With this prescription, the 
ratio cT(q, — q)/|q — q| evaluates to 1/2. 

Symmetrizing over ki and k2, and omitting ultraviolet power laws, we find 



(50(ki)50(k2)). 3 (27r)35(ki + k2)^^^ f-hHnk + ^k' + ---] (87) 
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Two derivatives. We have now exhausted all diagrams with only a single derivative. 
The next set of diagrams all involve two derivatives and break naturally into three sets: 
the first class includes all diagrams with the derivatives applied to both external legs 
of the two-point function; the second set includes all diagrams where one derivative is 
applied to an external leg while the other applies to the internal propagator; and the 
third set includes all diagrams with both derivatives applied to the internal propagator: 



(50(ki)50(k2)). D - Q . + + Q (88) 

Consider first the set of diagrams with both derivatives on the external legs. We 
obtain 

SL- : nD-^(l + gV)Q, (89) 
where Q is the quantity 

^ z(-ki,q) ■ z(-k2,q) 3 a(-ki, q)a(-k2, -q) _ q- z(-k2,q) 
^~ |q-ki|4|q + k2|2 ^4 |q-ki|2|q + k2p Iq-ks^ 

+ [ki^k2]. (90) 

Unlike the previous examples, none of the ratios which appear in Q are ill-defined. 
However, this result can still be simplified using the symmetry properties of z and a. 
In particular, we observe that cr is a quadratic form, and therefore 

(7(— a, — b) = (j(a, b) and z(— a, — b) = z(b, a). (91) 

These identities can be used together with the obvious antisymmetry of z [i.e., z(a, b) = 
— z(b,a)]. After performing the symmetrization over ki and k2, Q can be reduced to 
the simpler form 

z(q.k)2 3a(k,q)2 q-z(k,q) , , 

In this expression, k can be taken to be either ki or k2; after integration, the result 
depends only on the magnitude k and not its orientation, and we obtain 



(50(ki)50(k2))* D (27r)3(5(ki + ks)-^ i—kHnk--k^ + ---]. (93) 



' —k-^ Ink P 

The set of diagrams with one derivative on an external leg and one derivative on 
the internal propagator are the most complicated. To evaluate them, we write 

_J1 - 



: \l^- — {l-ikri){l-iqr])R, (94) 
Iq 

where R can be expressed as 

z(q.k)^ crfq. k)o'(k, q) 3 k ■ z(k. q) q-z(q,k) , , 

|q + k|6 |q + k|4 2 |q + k|4 Iq + kj'^ ^ ' 

in which we have used the symmetry properties of z and o", and the same convention 
that k may be chosen as either ki or k2 applies. In arriving at this expression for Q, we 
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have discarded a number of contributions of the form X- {z(— k, k) +z(k, — k)} for some 
vector X, which may itself require regularization. However, no matter how we choose 
to regularize the bracket {■■■}, the antisymmetry of z guarantees that it sums to zero, 
and therefore that such contributions cancel out of the observable expectation value. 
After integration, one obtains 

(50(ki)(50(k2))* D (27r)^5(ki + k2)^^^ (-kHnk+lk' + --] . (96) 



7r2 V 6 

The final class of diagrams of this type involve both derivatives applied to the 
internal propagator, 

which does not require regularization. After integration, one obtains 

(50(kO50(k2)). D (27r)=^5(ki + k^)^^ (i^^'l^^ - + ■ ■ ■) • (9^) 

Three derivatives. The only remaining class of diagrams are those containing three 
derivatives at the vertex. These diagrams break into two groups: those in which one 
end of the internal propagator is free of a derivative, and those in which an external leg 
is free of a derivative: 

(50(ki)50(k2)). D + (99) 

Both types give rise to comparatively simple expressions. For the first we obtain 

-J^: nD|^^(l-igr,)(4^i^ + l); (100) 
after integration this class of diagrams give contributions totalling 

P*{kf ( 1 3 1 3 



(50(kO50(k2)), D (27r)^5(ki + k2)-^ [--k^Tik^ -k' + ---j. (101) 
On the other hand, for the second type of diagram we obtain 

-J^: n2^,"(l-i*,)(4^ + l). (102) 

After integration, we find 

(50(ki)50(k2)). D (27r)35(ki + k2)^^^ (^kHTik- 1a;3 + . . .y (103) 
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5.3. Infra-red behaviour 

Having obtained the relevant contributions to the two-point function, given by Eqs. ( l85l) . 
(IHTj) . ( 1931) . (1961) . (1981) . ( TTOTll and ( ITOHD . we may collect these quantities to obtain an 
estimate of the total loop correction. It can be written 



where /3 is an unknown renormalization-scheme dependent quantity left over after 
cancellation of the ultraviolet divergences; it implicitly contains the (ultraviolet) scale 
which makes In k dimensionless. The coefficient of the logarithm, however, is scheme- 
independent [651 El] • 

6. Discussion 

In this paper, I have computed estimates for the leading radiative corrections to 
the one- and two-point expectation values of the inflaton field perturbation during 
a phase of single-field, slow-roll infiation. After suitable ultraviolet renormalization, 
the loop correction to the one-point function was found to be given by an arbitrary 
renormalization-scheme dependent constant which can be absorbed into the background 
value of (j). On the other hand, the loop correction to the two point function yielded a 
correction to the power spectrum of the form 



Although the amplitude of the 50 power spectrum itself is not observable, the 
amplitude of ( is accurately known to be of order 10^^°. At tree-level in a single- 
field model, the two are approximately related via ~ V^/e, where e is the slow- roll 
parameter introduced in Eq. ([5]). Since e is expected to be of order 10~^ or less, we 
can conservatively suppose that < 10~^°. The loop correction given by Eq. f llOSp is 
therefore extremely small provided that the logarithm is not too large. 

This does not yet allow us to conclude that loop corrections are too small to be 
observable in the CMB, because it is the loop corrections in ( rather than the 6(j) 
themselves which are accessible to experiment. Therefore, the prediction fllOSp must 
be translated into a prediction for 7?^'°°? before a final determination concerning the 
magnitude of loop corrections can be made. This calculation will be presented elsewhere 
[77] . However, it is already clear from Eq. fllOSp that quantum effects do not greatly 
disturb the magnitude of the fluctuations imprinted in 6(j) as successive k-modes pass 
outside the horizon. It is only the accumulation of fluctuations on superhorizon scales, 
where the flelds are in an effectively classical regime, which can give rise to a large loop 
correction. 

Eq. (llOSp is consistent with previous estimates which have been made in the 
literature. In particular, Weinberg has estimated a correction to from matter loops 
in a multi-fleld theory [65] which has the same functional form as fllOSp . Sloth [16| W7\ 




(104) 




(105) 
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has given a similar estimate, based on the same action given in Eqs. fl25p - fll7p . but 
evaluated several tens of e-folds after horizon crossing when large infrared divergences 
can compensate for a suppression in powers of slow-roll parameters; in this limit a 
different set of terms extracted from Eq. (fTTIl dominate the loop correction. This loop 
correction is proportional to (50^) ~ V* In(fc) for a flat spectrum, which reproduces the 
logarithmic /c-dependence described by (11051) . 
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